Eigenvalues of the Stewart— Lyth equation for inflation with a blue spectrum 
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By using the rather stringent nonlinear second order slow-roll approximation, we reconsider the 
nonlinear second order Abel equation of Stewart and Lyth. We determine a new blue eigenvalue 
spectrum. Some of the discrete values of the spectral index n 3 have consistent fits to the cumulative 
COBE data as well as to recent ground-base CMB experiments. 



I. INTRODUCTION 



For a wide range of inflationary models M, a single scalar field, the inflaton, is presumed to be rolling in some 
underlying potential U(<p). This standard scenario is generically referred to as chaotic inflation. A variance of it is 
the recently proposed quintessence. 

Our rather ambitious aim of reconstruction is to employ observational data to deduce the complete functional 
form {/(</>) of the inflaton self-interaction potential over the range corresponding to large scale structures, allowing a 
so-called graceful exit to the Friedmann cosmos O . The generation of density perturbations and gravitational waves 
has been extensively investigated. The usual strategy is an expansion in the deviation from scale invariance, formally 
expressed as the slow-roll expansion |J. 

In general, exact inflationary solutions, after an elegant coordinate transformation [Q, depend on the Hubble 
expansion rate H as a new "inverse time" , and the regime of inflationary potentials allowing a graceful exit has 
already been classified |^|. In our more phenomenological approach, the inflationary dynamics is not prescribed by 
one's theoretical prejudice. On the contrary, in this solvable framework, the 'graceful exit function' g{H), which 
determines the inflaton potential U{<p) and the exact Friedmann type solution, is reconstructed in order to fit the data. 
Even more, the transparent description of inflationary phase within the H— formalism was the foundation for deriving 
a new mechanism of inflation, called assisted inflation Q . 

Recent astronomical observations by COBE of the cosmic microwave background (CMB) confirm that the Universe 
expands rather homogeneously on the large scale. From cumulative four years CMB observations, the spectral index is 
now measured by COBE as n s = 1.2 ±0.3, including the quadrupole (/ = 2) anisotropy 0. This experiment has been 
complemented by ground-based CMB experiments. The preliminary data jq] from the Cambridge Cosmic Anisotropy 
Telescope (CAT) e.g., are consistent with COBE, albeit a slightly higher spectral index of n s = 1.3 ± 0.4. This 
allows H to constrain cosmological parameters, such as matter density fi m = 0.32 ± 0.08 and the Hubble constant to 
Hq = 47 ± 6 km/(s Mpc). A more recent analysis j[o| of the CMB anisotropy strongly constrain the spatial curvature 
of the Universe to near zero, i.e. k = 0. All data imply that a total density f2 = 1 (or k = 0) is 2 x 10 7 more probable 
then f2 = 0.4. Nowaday's best-fit to all CMB data is a Hubble constant of Hq = 65 km/(s Mpc) and requires dark 
energy with = 0.69 and a spectral index n s = 1.12. 

Stewart and Lyth [TT| proceeded from the exact power-law inflation in order to analytically compute the second 
order slow-roll correction to the standard formula jL2|. These analytical results are at present the most accurate 
equations available. This remarkable accuracy has recently been confirmed Jl| by numerical perturbations of the 
exact analytical result for the power-law inflation. 
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We apply the H -formalism to this accurate second order perturbation formalism, and transform the nonlinear 
equation of Stewart and Lyth into an Abel equation. In extension of our earlier work [|l4| , we determine here a discrete 
eigenvalue spectrum of this nonlinear equation. In particular, our new blue discrete eigenvalue spectrum has the 
scale invariant Harrison-Zel'dovich spectrum with index n s — 1 as a limiting point. Thus some specific eigenvalues 
1 < n B < 1-5 of the family of solutions can be related rather accurately to existing as well as future observations. 



II. GENERAL METRIC OF A SPATIALLY FLAT INFLATIONARY UNIVERSE 

In the following we reconsider the Lagrangian density of a rather general class of inflationary models in the Einstein 
frame: 



C E = ^y/[g~\(R + K [g^(d^)(d^)-2U^)}) . (2.1) 
in 



Here R is the curvature scalar, </> the scalar field, t/(</>) the self-interaction potential and k — 8ttG /c 4 the gravitational 
coupling constant. We use natural units with c = h = 1 and our signature for the metric is (+1, —1, —1, —1) as it is 
common in particle physics. 

For the flat (k = or O = 1) Robertson- Walker metr ic fa vored by the inflationary paradigm and by recent obser- 



vations, the evolution of the generic inflationary model (2.1) is determined by the autonomous first order equations 



H = g(<l>,H), </> = ±y/2/Ky/-g(<l>,H), (2.2) 

where g(<f),H) := nU((f>) — 3H 2 . De Sitter inflation with H = is a singular, but well-known subcase, cf. [pjjfl . For 
g(4>,H) =/= 0, however, we can simply reparametrize the inflationary potential U((f>) — ► U = U(H) := U(<j)(t(H))) 
in terms of the Hubble expansion rate H := a(t)/a(t) as the new "inverse time" coordinate. For this non-singular 
coordinate transformation, the term g((f), H) reduces to the graceful exit function g(H) = kU(H) — 3H 2 < and the 
following general metric and scalar field solution of (2.2) emerges ||,[|: 



, o dH 2 o f f HdH\ r ., t ) 

ds = ~T77T\ - «o exp 2 ' > Fj ~ 2 1 ~ 2 1 Aa2 1 ~ : - 2 



g 2 (H) r V J g(H) 



j [dr 2 +r 2 (d6 2 + sin 2 6chp 2 )] , (2.3) 



Once we know g(H) from some experimental input, we can (formally) invert ( p.4|) in order to recover from the chain 
of substitutions 

[g(H) + 3H 2 ]/k = U(H) = U(H o 4>- x {4>)) = U(<f>) (2.5) 
the inflaton potential in a simple and elegant manner. 



III. SECOND-ORDER SLOW-ROLL APPROXIMATION IN CANONICAL ABEL FORM 



For a reconstruction of the potential U ((f)) under the "umbrella" of chaotic inflation, we adopt here the equations 
of Stewart and Lyth fllf for the second order slow-roll approximation. In terms of the energy density e = —g(H)/H 2 , 
this nonlinear second order differential equation involving the scalar spectral index n s reads p4|Jlq| 



where -=d/d\nH 2 , C 



2Ce e -(2C + 3)ee - e +e 2 + 
-2 + ln2 + 7~ -0.73, and 

A := (n s - l)/2 







(3.1) 



(3.2) 



denotes the deviation from the scale invariant Harrison-Zel'dovich spectrum. The corresponding second order equation 
for the spectral index n g of tensor perturbations turns out to be 



n g = — 2e 



1 + e - 2(1 + C) e 



1 . 



(3.3) 
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where the approximation is valid in first order. Note that for the constant solution e — Aq, where 



A 2 + A + A = 0, => A = i (-1 T V3 - 2n s ) , (3.4) 
we recover from the second order equation exactly the first order consistency relation n g = n s — 1. Moreover, 



we can infer already from (3.4) that a real eigenvalue spectrum is constrained by n s < 1.5. 



By introducing the flow of the energy density e via p(e) :=e with e= pdp/de, the second order equation (3.1) can 
be brought to the following form: 

- (P + «e)p + /3(e 2 + e + A) = , (3.5) 

de 

where a = 1 + 3/(2C) and /3 = 1/(2C). For invertible energy flow, it can be transformed via u := 1/p into the 
canonical form of the Abel equation of the first kind, cf. |17| , 

nil 

2Ce— + u 2 [1 + (2C + 3)e] - u 3 [e 2 + e + A] = . (3.6) 
Garcia et al. |l8|] have demonstrated that, within the class of finite polynomials p = X)i=o ^ e *' ^ ne um< 4 ue solution 



of (3.5) is given by a polynomial p = A + b\t of /irst degree, where b\ = (a± y/a 2 — 4/3)/2 necessarily satisfies the 
additional constraint 6i (1 — 6i ) = A due to the nonlinearity of the Abel equation. Mathematically, this transformation 
improves the order of the polynomial by one in comparison with (3.1), where the solution is just a polynomial of zero 
degree. This first order polynomial corresponds to e = A a + e y 1+A ° : where A is again given by ( |3.4[ ). Since A is 
determined due to the additional constraint mentioned above, cf. the spectral index n s is now completely fixed 
by the Euler constant 7 as 

2C + 9 T 3V4C* 2 + AC + 9 
= 1 ic3 

= 1.49575 (-6.57797) . (3.7) 

In order to circumvent this 'no-go' theorem for polynomial solutions, we shall introduce a different transformation 
which leads us to new solutions in form of an infinite series. 



A. New class of solutions with a blue spectrum 

If we assume that e = e(y) with y := H 2 , Eq. ( |3.l[) can be rearranged to 

2Cey 2 e" - (3e + l)ye' + e 2 + e + A = , (3.8) 

where '=d/dy. We try the ansatz 

00 

e = ]T^y\ (3.9) 

8=0 

and find that the odd powers have to vanish. This can be traced back to the fact that ( |3.8[ ) is invariant under 
reflections y — > — y only for even functions e(— y) = +e(y). For even powers, the following Taylor expansion can be 
recognized: 

00 

e = J2 A 2mjy 2mj , (3.10) 
j=o 

where m is a positive integer, which labels the first non-constant term in our expansion and thus distinguishes the 
different subclasses of solutions. Hence, the series for m = 1 starts as e = Aq + A2J/ 2 + + . . ., for the case m = 2 
as e = + A4J/ 4 + Ag?/ 8 + ■ . ., a nd so on. In all cases, the coefficient Aq of zeroth order (y°) is again determined 
algebraically by a A via Eq. ( |3.4[ ); or, as we shall show below, Aq is determined by the next higher order, and the 
zeroth order establishes n s . 

In order 2to, we find the following discriminating relation 
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{[2m(2m - 1)C - 3m + 1]A - m + 1/2} A 2m = , (3.11) 

due to the nonlinearity of the Abel type equation. This recursion relation implies Aim = (and e = Aq = constant) 
or the new additional constraint 

A ° = 2[2m(2m-"l)C' 1 -3m+l] < ° ' (3 ' 12) 



Together with (3.4) this implies that the spectral index of our new class of solutions can only adopt the following 



discrete eigenvalues 



(2m-l)[4m-l-4m(2m-l)C] 

2[2to(2to- l)C-3m + l] 2 ' 1 ' ' 



For solutions with large to, we find asymptotically 

^lir "-^"db' (3 - 14) 

such that the scale invariant Harrison-Zel'dovich solution with index n s = 1 turns out to be limiting point. Let us 
stress that our new discrete spectrum approaches it from the blue side, which previously was considered rather difficult 
to achieve. In order to correlate this with observational restrictions, let us display the highest eigenvalues for n s . For 
m = 1, we have (ACAq — AAq — 1)A 2 = and, because of A 2 ^ 0, 

A ° = 4(C - 1) = ~°' 1445 ' (3 - 15) 



Eq. (3.4) discovers the value of the scalar spectral index to be 

n s = 1.247 . (3.16) 

The integration constant A 2 remains arbitrary. The fourth order establishes the constant A4 as an algebraic function 
of A , A 2 and similar for higher orders. 
For to = 2, the fourth order constraint is 

^=24^10 = - °- 109 ' (3 ' 17) 
while A4 now plays the role of a free constant. Then, zeroth order vanishes if 

n s = 1.194 , (3.18) 

which is very close to the experimental value n s = 1.2 ± 0.3 from COBE. Furthermore, it is revealed that As = 
Aq), A\ 2 — f 2 (As,A4, Aq), and so on (where fi,f 2 are some algebraic relations). The higher order coefficients 
of the series can be derived recursively via Mathematica, similarly as in the case of the Bartnik-McKinnon solution, 
cf. |3. 

IV. RECONSTRUCTION FOR THE FULL EIGENVALUE SPECTRUM 

In cosmological applications, the energy density e in our second order reconstruction has only a limited range. For 
the decelaration parameter q(t) := — a a/ a 2 = e — 1 to become negative, as necessary for inflation, our solutions are 
constrained by < e < 1. In the case of the constant solution e — Aq, merely the spectral index is restricted to the 



continuous range —3 < n s < 1 . Since this corresponds to g(h) = —AqH 2 , a simple inversion of the integral (2.4) and 
subsequent insertion into ( |2.5|) leads to the potential 

U{cj>) = ^-^C 2 exp (±^A~ 0) , (4.1) 

of power-law inflation, where C is an integration constant. 
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Similarly, for the ex trem e blue spectrum (3.7) the graceful exit function reads g{H) = —H 2 \Aq + €qH 
yielding via (2.4) and (2.5) the reconstructed potential 



2(1+A 



A(i + tan 2 (P0)) 



1/(1+A ) 



A o tan 2 (P0)) , 



(4.2) 



where P = ±-v/— kAq/2 (1 + A ) is a constant. According to its shape, cf. |19|, it interpolates between the power-law 
inflation in the limit eo ~ > and intermediate inflation. Since n s ~ 1.5, the parameters are approximately Aq ~ — 1/2 
and P ~ y^/4- Its form is similar to the exact solution of Easther |H| except that the latter yields the too large 
value n s = 3 for the spectral index. 

For the new bl ue spe ctrum (3.13) we have reconstructed the potential U (<ft) numerically in Fig. [I] and || for the first 
two indices ( [3.16 3.18| ), see the Appendix for details. The form of both new potentials resemble the ones found in 
[^l|,0. The corresponding forms of kII(H 2 ) are shown in Fig. || and || Inflation occurs if the scalar field is between 
the two lines 2P 2 and 3P 2 , cf. {§. 

The range of 4> for which these potentia ls a rc valid depends also on the (so-far missing) information on the spectral 



index n g of tensor perturbations, cf. Eq. (3.3). For our solutions for to = 1 and m = 2 we show the cf> dependence of 

Aq, we expect the first 
-1 G (-0.4,0.5] @. In 



m Fig. | and §. Since all our new solutions are 'deformations' of the constant solution e 
order consistency relation of power-law inflation to hold to some approximation, i.e. n g ~ n 
second order, we recognize from Fig. [5] and |] that for our solutions n g is always negative. 

Of particular interest j2^] is also the squared ratio of the two further observables, the amplitude A g of gravitational 
waves, versus its counter part A s , the amplitude of primordial scalar density perturbations: 



(A 2 g /A 2 s 



l-2Ce 



~ (V3 - 2n 8 - 1) 



(4.3) 



where the approximation is valid in first order. We display the numerical results in Fig. and [3] as well. Form first 
order approximations, we expect the range of this ratio within ~ [0, 0.2). For several well-known inflationary models 
(power-law and polynomial chaotic, e.g.) the predicted ratio can be found in Table 2 of p2| . Recently, assisted 
inflation was able to produce a scalar spectral index closer to scale-invariance (but for n s < 1) because of its large 
number of non- inflationary real scalar fields present in this theory || . 

In this paper, we followed the basic idea of reconstructing the inflaton potential in second order, given in Q. But, 
instead of producing the potential just for a few values <f>, we could construct the complete functional form of the 
potential for a series of constant spectral indices. As it is well-known, a full reconstruction can only be achieved if 
informations on the gravitational parts are included; scalar perturbations determine the potential up to an unknown 
constant. Since the differential equations are nonlinear, different constants yield different potentials but for the 
identical scalar spectrum. Any piece of knowledge concerning the tensor perturbations breaks this degeneracy; so we 
use our solutions to calculate both the tensor index n g and the ratio of the amplitudes A 2 /A 2 necessary in second 
order to define the inflaton potential uniquely. We recognize that both parameters are not constant in contrast to 
the scalar spectral index. We have checked numerically several consistency equations in second order and found good 
confirmation within the inflationary region. Furthermore, we analyzed the behavior of the slow-roll parameters e and 
j] = —dg/dH 2 , cf. the definition in At the beginning of inflation, we derive that rj > 2e which is valid for the 

(to = l)-solution up to e ~ 0.32 and for the (m = 2)-solution up to e ~ 0.28. This feature resembles that of hybrid 
inflation p3j , cf. Table 1 of , where the influence of gravitational waves is negligible in accordance with the small 
values of both n g and A 2 /A 2 as can be recogn ized within Fig. || and|| close to 4> ~ 1. 

Summarizing, the nonlinear Abel equation (3.1) permits exact solutions for the discrete or continuous eigenvalues 
of Table H. All positive discrete values are marginally consistent with the n s — 1.2 ± 0.3 from COBE as well as recent 
ground-based CMB experiments, our highest value for the blue spectrum is just below the n s < 1.5 constraint in 
order to suppress distortions of the microwave spectrum and the formation of primordial black holes during the phase 
of reheating, cf . B . 



V. STRING OR M THEORY AXION AS QUINTESSENCE 

Recently, quintessence in the form of slow-rolling scalar field Q has been proposed as an alternative form of dark 
energy density with negative pressure pi] , and Refs. therein. The true minimum of the quintessence potential is 
presumed to vanish, i.e. (J7q) m i n = 0, however the present value of Q is displaced from the true minimum, providing 

U Q - (3 x 10- 3 eV) 4 , (5.1) 



5 



and a negative pressure with the equation of state 



gQ = Q 2 -2£/ Q < 
PQ Q 2 + 2U Q 



while obeying the usual equations of motion. 



Our potential (4.2) for the blue spectrum (3.7), mimics to some extend the required properties of quintessence, 
cf- g]. 

Given the assumption of (i/Q)min = 0, one can in principle compute the quintessence potential Uq once the particle 
physics model for Q is given. However, the required size of Uq is extremely small compared to any of the mass scales of 
particle physics, and thus an utmost question is how such a small Uq can arise from realistic particle physics models. 
In fact, the most natural candidates for a light scalar field whose typical range of variation is of order Mpi anc k/y§7r 
form string or M -theory moduli multiplets describing the (approximately) degenerate string or M-theory vacua [p6| . 
It is thus quite tempting to look at the possibility that Q corresponds to a certain combination of the string or 
A/-theory moduli superfields. 

Current ground-based CMB experiments are already discriminating between different inflationary models and 
are preliminarily constraining the cosmological parameters. In future, the satellite-based microwave background 
anisotropy measurements from PLANCK Surveyor will constrain J27[ them even more, likely to less then 1 %. 
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APPENDIX: DERIVATION OF THE NEW POTENTIALS 

Within the if-formalism, real scalar fields can only be determined if kU(H) < 3H 2 [J|||. Furthermore, inflation 
just occurs within the area 2H 2 < kU(H) < 3H 2 . We use this knowledge to calculate numerically the complete 
potential in real <p space. For instance, for m = 1 we define a new function f(y) by 

e =: A + f(y)y 2 (Al) 

so that its corresponding differential equation has the initial values /(0) = A2 and /'(0) = 0. The regime of real 
solutions can be derived from the differential equation for the scalar field 

0' = tJI^= ■ (A2) 
V K ^/4ey 2 

Because Aq < we shall find an area around y — where <f> is imaginary, hence, nU(H) > 3H 2 . The point defined 
by y = -J— Ao/f(y), corresponding to kU(H) = 3H 2 , establishes then the initial value f(y) (if y is chosen) in real <fi 
space. f(y) is still connected with the initial value A2 in complex <j) space. The potential can be obtained immediately 
from the differential equation 

kU' = (3 - A Q ) - y 3 f - 3y 2 f . (A3) 

In order to derive Fig. [|, we used A2 = 0.01 and found y « 0.008026 which is linked to f(y) « 2250. Furthermore, 
we applied the initial values f'(y) = —1 and 4>(y) = 1. The initial value for the potential is, of course, nU{y) = 3y. 
By this, the point <f> = 1 in Fig. [j] is correlated to k,U(H) = 3H 2 (no potential does exist for (/)>!); inflation 
(kU(H) > 2H 2 ) occurs till almost the maximum of the potential; for smaller values of <fi kinetic energy dominates 
and the scalar field becomes constant. We have investigated a broad interval of f(y) (between 0.1 and 10 6 ) and (f>(y) 
values but always found this form of the potential. For A2 < 0, the function f(y) is always negative, hence, </> always 
imaginary. For m = 2, we put f[y) = 2250 so that y m 0.083427; f(y) = -1 and tf>(y) = 1. 
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FIG. 1. Inflationary potential U(4>) with n s = 1.247 (m = 1). Further explanation in the Appendix. 
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FIG. 2. Inflationary potential U(<j>) with n s = 1.194 (m = 2). Further explanation in the Appendix. 
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FIG. 3. Inflationary potential nU(y = H 2 ) for in = 1 (potential function until y = 2.4) and m = 2 (until y = 17). The two 
lines correspond to kU = 2y and kU = 3y, the boundaries of the inflationary phase. Further explanation in the Appendix. 
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FIG. 4. Enlarging the surrounding of the origin of Fig. |3j. It is recognizable that both potentials start at the line kU = 3y, 
below the boundary where real scalar fields are present. For the (m = l)-solution, the initial value y is approximately 0.008026 
(for the values given in the Appendix), and for m = 2, y as 0.083427. 
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FIG. 5. The ratio of the tensor and the scalar spectrum A 2 g /A? s and the tensor spectral index n 9 for the (m = l)-solution 
within inflationary region. 




FIG. 6. The ratio of the tensor and the scalar spectrum A g /A s and the tensor spectral index n g for the (m = 2)-solution 
within inflationary region. 



Eigenvalue n s 



TABLE I. Eigenvalues of the spectral index 



Spectrum 



Potential U{<j}) 



(-6.57797) 

-3 < n s < 1 



1.106 
1.125 
1.153 
1.194 
1.247 
1.49575 



(unrealistic) 
continuous 

discrete 
(m = 5) 
All CMB data 
(m = 3) 
COBE data 
(m = l) 



~ tan^v^yAl) 
oc exp(±v / 2K^4o <j) 



cf. Fig. 2 
cf. Fig. T 

tan 2 (yK(?!)/4) 
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